Quantum Hall effect in narrow graphene ribbons 
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The edge states in the integer quantum Hall effect are known to be significantly affected by electrostatic 
interactions leading to the formation of compressible and incompressible strips at the boundaries of Hall bars. 
We show here, in a combined experimental and theoretical analysis, that this does not hold for the quantum Hall 
effect in narrow graphene ribbons. In our graphene Hall bar, which is only 60 nm wide, we observe the quantum 
Hall effect up to Landau level index k = 2 and show within a zero free-parameter model that the spatial extent 
of the compressible and incompressible strips is of a similar magnitude as the magnetic length. We conclude 
that in narrow graphene ribbons the single-particle picture is a more appropriate description of the quantum Hall 
effect and that electrostatic effects are of minor importance. 

PACS numbers: 73.63.-b, 81.05.ue, 73.43.-f 
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I. INTRODUCTION 

The discovery of graphene is intimately related to the first 
measurements of the integer quantum Hall effect in this mate- 
rial system^* 2 - It was observed as early as 2005 that the inte- 
ger quantum Hall effect shows an anomalous sequence of Hall 
plateaus in the Hall conductivity o~ xy = 4(fe + 1/2) e 2 /h due 
to the graphene-specific Berry phase that finds its origin in the 
bipartiteness of the honeycomb lattice. This is a crucial dif- 
ference from quantum Hall effect measurements in more es- 
tablished material systems such as two-dimensional electron 
gases (2DEGs) in semiconductor heterostructures. 

A second distinction between the quantum Hall effect in 
graphene and that in ordinary 2DEGs is suggested by two ter- 
minal conductance studies on narrow ribbons,^ where con- 
ductance quantization was interpreted in quantum Hall terms, 
even though the ribbons are narrower than the typical separa- 
tion between edge states in a 2DEG. To confirm the validity of 
this interpretation, we have fabricated a 60-nm-wide graphene 
Hall bar, allowing for proper four terminal measurement of 
R xy and confirming that quantum Hall edge states do survive 
even in such narrow devices in graphene. 

In order to reconcile these observations with the conven- 
tional picture of edge-state separation, we need to recall the 
physics of edge states in ordinary 2DEGs. The first detailed 
description of edge-state transport in the quantum Hall regime 
was provided by HalperinS, using a noninteracting electron 
picture. While this model is useful in describing much of the 
relevant physics, it was subsequently found that such a pic- 
ture does not in general properly account for the position of 
the edge states. A more exact description emerges when con- 
sidering the influence of electrostatic interactions between the 
electrons that form the quantum Hall state and the surround- 
ing electrons in the substrate and the nearby metallic gates. In 
ordinary 2DEGs, it is well established that these electrostatic 
interactions lead to the formation of compressible and incom- 
pressible strips^ Indeed, the effect of electrostatic interac- 
tions is of such importance that, in many cases, the single- 
particle picture fails to provide even a qualitative description 
of experiment. The situation is different in narrow graphene 




FIG. 1 . (Color online) SEM image of a plasma patterned six-terminal 
graphene Hall bar with w = 60 nm and I — 250 nm. 



structures, where, as our zero-free-parameter model shows, 
the relevant length scale for the single-particle picture — the 
magnetic length Ib = \J h/eB — is similar to the width of 
the compressible and incompressible strips. 



n. EXPERIMENTAL PART 

A scanning electron micrograph (SEM) of our sample is 
shown in Fig. [I] It is made from single-layer graphene flakes 
exfoliated from an HOPG block onto a highly p-doped Si sub- 
strate having a 285 nm thick Si02 cap layer. After contacting 
Ti/Au electrodes, these flakes are etched in an Ar/02 plasma 
to obtain narrow Hall bars with six terminals. To reduce resid- 
ual contaminations, the etched structure is annealed at 200 °C 
in forming gas for about 2 h. The sample is studied at the base 
temperature of a dilution refrigerator using standard low- noise 
ac measurement techniques. 

Preliminary characterization of the sample is done at 4.2 K, 
where, as typical for ribbons of width below 100 nm, the con- 
ductance G of the sample is strongly suppressed in the vicinity 
of the charge neutrality point, which is at V® ~ 7.5 V. At high 
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FIG. 2. (Color online) Hall conductivity a xy as a function of gate 
voltage V g and Landau level filling factor v — h/(e 2 B]C g (V g — 
V g u ) at B = 11 T and 30 mK. (Figure adapted from Ref. @.) 

electron and hole concentrations, the field-effect mobility of 
the charge carriers is fx = 2a 'c 1V~ ~ 1500cm 2 /Vs, where 
the effective gate capacitance C g « 210/iF/m 2 is determined 
from Hall measurements at various gate voltages V g . This 
corresponds to a diffusion constant of D ~ 0.01 m 2 /s and a 
mean free path of l m ~ 20 nm. Such small values are typ- 
ical for plasma-etched graphene nanoribbon devices and are 
a consequence of the great amount of edge disorder in these 
systems^ 

We recently observed the quantum Hall effect in very nar- 
row graphene wires.— A representative example of such a mea- 
surement is shown in Fig. [2] where we plot the measured Hall 
conductivity a xy as a function of the back-gate voltage V g and 
filling factor v = ±4(/c + 1/2) at B = 11 T. The occurrence 
of conductance plateaus at k = 0, 1, and 2 is typical for Dirac 
fermions and clearly confirms that our sample is mono layer 
graphene. In general, a xy plateaus are well developed if the 
energy spacing between the corresponding Landau levels is 
much larger than their broadening T. In graphene the former, 
given by Eq. (Q, decreases with increasing filling (i.e., the in- 
dex fc), whereas the broadening T = hvp/l m ~ 30meV is 
constant. This explains why only the k = 0, 1, and 2 plateaus 
are visible up to B = 11T, with the k = plateaus being 
most pronounced. We found this observation remarkable in 
view of previous experience with (Al,Ga)As-based structures. 
Therefore, we decided to investigate the electrostatics of the 
edge channels in more detail, which is presented in the next 
sections. 



IE. THEORETICAL PART 

From a theoretical point of view, the essence of our device 
is an etched graphene ribbon separated from a metallic back 
gate by an insulating Si02 layer. We note that the thickness, of 
the gate insulator is a significant parameter because, for very 
thin thicknesses tunneling becomes as important as electro- 



static effects. Our geometry is therefore fundamentally differ- 
ent from the one used in recent STM studies of the spacial ex- 
tent of edge channels in graphene resting directly on graphite, 
where the thickness of the gate insulator is basically (i.e., the 
distance between the graphene and the graphite substrate) J£ 
We now examine the consequences of our geometry on 
the accepted models by Chklovskii, Shklovskii, and Glaz- 
man (CSG} 6 as well as Chklovskii, Matveev, and Shklovskii 
(CMS)jZi In the electrostatic theories of CSG and CMS, a 
2DEG is laterally confined by electric fields originating from 
gate electrodes which are in plane with the 2DEG. In the ab- 
sence of magnetic fields, the electrostatic solution of such 
a two-dimensional problem was provided by Larkin and 
Shikin.'i They find a charge density distribution which re- 
duces from its bulk value to at the boundary of the electron 
gas. Although the confinement of mobile electrons in these 
models is generated by gates, the authors of Refs. |6) and 7 
argue that their results also apply to etched structures of con- 
ventional heterostructures. This is justified by surface states 
which lead to Fermi level pinning and negative charges accu- 
mulate at the edges. Inner electrons are repelled, yielding the 
same depletion effect as split gates. However, as shown by 
Silvestrov and Efetov (SE),— the situation is substantially dif- 
ferent in the case of a graphene nanoribbon. Since the edges of 
this two-dimensional material end abruptly, they act as hard- 
instead of soft-wall boundaries. Assuming a constant electro- 
static potential in a ribbon positioned in the xz plane, such 
that in the transverse x direction no electric force acts on its 
excess electrons, SE derived an expression for the charge den- 
sity. If the ribbon width 2a r is smaller than the gate dielectric 
thickness b, the electron density across the ribbon is, to a good 
approximation, given byi^ 




(l) 

where Co = 7^°° i -r— I tt^ ] t^t is a numerical 

constant 1 ^ that, for a dielectric constant e ~ 3.9 as in the case 
of SiC>2, has a value of 0.175. The charge per unit length a is 
linearly related to the back-gate potential V 9 ,— 

TT§) Inn- 

The inverse-square-root edge singularity in Eq. (H) is in 
stark contrast to the soft density profiles (p ~ ^fx) predicted 
at the edges of conventional heterostructures . 6 ' 7 ' 1 1 In the fol- 
lowing, we analyze the impact of this diverging density distri- 
bution on the formation of compressible and incompressible 
strips in graphene. 

The formulas above are valid for any metallic strip with 
a hard-wall boundary. The graphene-specific physics comes 
into play with quantum effects. Evidently, the Pauli principle 
prevents the excess electrons from all being in the same state. 
When the Fermi wavelength varies slowly on the considered 
length scale (here the width of the ribbon), one can divide the 
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FIG. 3. (Color online) Schematic of the compressible and incompressible strips in graphene, showing the electron density oc n(x) and the 
Landau levels with their local fillings. In (a) the ribbon is in a C state with a compressible strip at its center whereas in (b) it is in an I state with 
an incompressible strip in the middle. The regions near the edges of the ribbon where the Landau levels increase very sharply due to hard-wall 
confinement are not shown. In these regions of size ~ Ib away from the edges, the electrostatic analysis is no longer applicable*^ 



system into many elements of sizes where equilibrium ther- 
modynamics is still applicable. By using this Thomas-Fermi 
approximation, SE defined a local Fermi momentum for each 
element and derived a relation between the inhomogeneous 
charge density given by Eq. ([U and the potential felt by an 
individual electron. At zero temperature and in the limit of 
rather narrow ribbons (a r <C b), the curly bracket in Eq. (Q~|i 
reduces to 1, and the potential has the form 



U(x) 



-Hvf 




(3) 



with vf the Fermi velocity in graphene. The derived expres- 
sion is valid at distances Sx > a ? V 3 (f) 2 ^ 3 from the strip 
edges, where the Thomas-Fermi approximation is applicable. 
This potential should be viewed as a quantum correction to the 
electrostatic potential and describes locally the energy differ- 
ence between the constant Fermi energy and the Dirac point. 

Knowing the electron distribution and the potential in a 
graphene ribbon, we can now adapt the CMS formalism to 
graphene. Making use of the semiclassical electrostatic pic- 
ture described above, we expect, for a finite magnetic field, 
the scenario sketched in Fig. [3] As in the case of the 2DEG 
in Ref.Q, the electron density ~ n(x) in our graphene ribbon 
is divided into alternating compressible and incompressible 
strips. The latter are characterized by filled Landau levels and 
describe local unscreened regions in which electrons have no 
possibility to reorder. This is due to the large energy gap they 
must overcome in order to find unoccupied states in one of the 
nearby empty levels. On the other hand, compressible strips 
are characterized by partially filled Landau levels residing at 
the Fermi level. The many unoccupied states allow for the 
electrons to rearrange and minimize their energy such that the 
electric field is locally screened. In general, two situations 
are possible. In Fig. [3ja), we show a C state that describes 
a central compressible strip in the fc-th Landau level which 
is partially filled near x — 0. Upon decreasing the magnetic 



field strength or increasing the density of excess charges, the 
fc-th Landau level gets filled and the ribbon is driven into an I 
state with an incompressible strip at the center [cf. Fig. 0b)]. 



IV. RESULTS 

The basic formalism used in Ref.0to determine the widths 
and positions of the compressible and incompressible strips in 
a gate confined 2DEG is also valid in graphene nanoribbons. 
Imposing the same requirements on the electrostatic boundary 
conditions as within the framework of the CMS model, and 
solving for our geometry leads to the following set of equa- 
tions: 



K0)-4(fc + l/2)]n L +n"(0)- 



(4) 



2 7ren"(0)a 
3(e+l) 



.U « J, 

rri (a 



+ a' 2 )E \yi - (a' /a) 2 - 
(a') 2 K [yi - (a'/a) 2 ] } = ~A*k > (5) 



TABLE I. Experimental parameters used in the modeling. 



Quantity 


Value 


Magnetic field B 


11T 


Insulator dielectric constant e 


3.9, Si0 2 


Thickness of insulator b 


0.3 /im 


Half-width of graphene ribbon a r 


30 nm 


Back-gate potential V g 


-20 to 35 V 
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FIG. 4. (Color online) Top: Two-terminal conductance plateau cor- 
responding to Landau level index k — 0, plotted over the back-gate 
voltage V g . The charge neutrality point is at Vf, and the dashed line 
indicates Go- Bottom: The position and width of the compressible 
and incompressible strips are plotted, which, in Fig. [3] correspond to 
the occupation level 4(k + 1/2) ul with k = 0. The dashed-dotted 
line shows a, the outer edge position; the solid line shows a , the in- 
ner edge position of the incompressible strips. In the upper graph the 
suppressed contribution of distant incompressible strips to the con- 
ductance at the I-C transition is incorporated, whereas in the lower 
graph it is neglected, resulting in a slight mismatch in width of the 
shaded areas. 



v(0)-v B (0) 



n"(0)(a 



An L 
(a -a') 2 



i>2 



1/(0) - A(k + 1/2)] . (6) 



Here, E and K are elliptic integrals, and the quantities a and 
a' determine the positions and widths of the compressible and 
incompressible strips next to x = 0. The energy difference, 



fJ-k 



hvp 
~Tb~ 



(7) 



is that of adjacent Landau levels (see Fig. [3) and depends on 
the index k, the Fermi velocity vp, and the magnetic length 
Ib = y/h/eB. The occupation number of these levels is 
defined by vb(0) = Ub{x = 0)/til, with ns(0) the elec- 
tron concentration at magnetic field strength B in the center 
of the ribbon, and jil = 1/2tt1 2 b the Landau level degener- 
acy. Furthermore, we define i/(0) = n(x = 0)/ul, with 
n(x) the electron concentration in the absence of a magnetic 
field, which can be extracted from p{x, y) = en(x) S(y) via 
Eq. A}. 



Taken together, Eqs. ©, (0, and © represent a com- 
plete system from which the quantities a, a', and ^b(0) can 
be derived. The latter is linearly related to the two-terminal 
conductance 7 via 



G=-i/ B (0). 



(8) 



In order to analyze the results provided by this model, we use 
the experimental parameters listed in Table [j] and solve for a, 
a', and vb- Since Eqs. (01, ©, and (O are derived under the 
assumption that n(x) ~ n(0) + ^n"(0) x 2 , we are limited 
to a parameter range where this approximation is valid. This 
is ensured when a and a' are small compared to a r which 
is the case at the transition from a C to an I state caused by 
increasing the doping or decreasing the magnetic field. In 
the upper graph in Fig. |4] we show the conductance plateau 
which corresponds to the Landau level with index k = 0. 
From the electrostatic point of view, incompressible strips 
next to x = 0, characterized by local regions with an addi- 
tional charge accumulation, are described by dipolar strips,^ 
Although the calculated charge density distribution, Eq. ([T), 
is not altered on a rough scale at any reasonable magnetic 
field, the formation of such strips leads to a small enhance- 
ment in charge concentration at x = 0. Thus, it is reasonable 
that the conductance, given by Eq. ([8]), lies above the dashed 
line, which indicates Go = T-f(0). From the lower graph in 
Fig. |U we can read off the positions and widths of the com- 
pressible and incompressible strips which correspond to the 
occupation level 4(fc + 1/2) Ul with k = 0. In the C state, 
before the formation of the conductance plateau, the incom- 
pressible strips next to the central compressible one have a 
width of ~ 9 nm ... 14 nm. This, however, contradicts a cru- 
cial assumption of the electrostatic model,— 



h^k < (a- a') . 



(9) 



For the quasi classical electrostatic treatment above to be jus- 
tified, the size of the electron wave function determined by the 
magnetic length should be much smaller than the length scales 
provided by the widths of the compressible and incompress- 
ible strips. Since Ib=iit ~ 8 nm, our experiment is in the 
crossover regime to the conventional single-particle picture of 
ballistic transport^ The latter is mandatory at high quantum 
numbers k where condition (0 is strongly violated. For the 
first several plateaus, however, the electrostatic model which 
results in Eqs. ®, (O, and (|6]l is still relevant, as Eq. © is 
approximately valid. Also note that the widths of incompress- 
ible strips located at xt scale as ~ 1/ y/n'(xk)-— Thus, for 
narrow ribbons with large n'(xk) and therefore small a — a', 
condition (0 is violated much earlier as for wide ribbons with 
expanded incompressible strips. Furthermore, since n'(xfe) 
increases as we approach the edges, the spatial extent of the 
edge states within the electrostatic picture is much narrower in 
our graphene ribbon compared to Ib ■ Hence, we show that the 
CSM picture of transport in the quantum Hall regime breaks 
down here and the simpler single-particle picture not only suf- 
fices but, indeed, is essential to correctly describe the physics. 
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V. COMPARISON OF THEORY AND EXPERIMENT 



Comparing the data on the k = plateau, corresponding 
to filling factor v = 2 in Fig. [2] with the theoretical result 
given in Fig. [4] we find discrepancies in width and location, 
both measured in units of the back-gate voltage. Our electro- 
static model leads to very narrow plateaus occurring at gate 
voltages which are too small by a factor of ~ 3 to 4. This dis- 
agreement is consistent with theory overestimating the gate 
capacitance, which is a commonly observed discrepancy in 
graphene nanoribbons;- The narrowness of the theoretically 
derived plateaus compared to their experimental widths can 
be explained by disorder, which is certainly present in the de- 
vice, but neglected in our model. As realized in Refs. |^ and 
7, disorder leads to localization in a compressible liquid a low 
small enough density, such that at the I to C state transition 
there is a range of voltage where Landau levels which locally 
(in the vicinity of x = 0) start to fill up cannot contribute to 
the conductance. 
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VI. CONCLUSIONS 

We have performed four terminal transport measurements 
on a 60-nm-wide Hall bar of graphene. Our results vali- 
date the interpretation that the Hall effect is responsible for 
the quantized conductance in narrow ribbons. A quantitative 
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implication for the spatial extent of edge states in such de- 
vices, confirms that, indeed, the formation of compressible 
and incompressible strips which dominate the physics of tra- 
ditional 2DEG devices is of minor importance here. Instead, 
the single-particle picture is more appropriate. 
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